MINIMUM VARIANCE UNBIASED AND MOMENT ESTIMATORS OF 
CARRIER FREQUENCY OFFSET IN MULTI-CARRIER SYSTEMS 

CROSS-REFERENCE TO RELATED APPLICATION 
This application is related to a co-pending application entitled "Globally 
Optimum Maximum Likelihood Estimation of Joint Carrier Frequency Offset and 
Symbol Timing Error," U.S. Ser. No. 09/496,890, filed on March 2, 2000, assigned to the 
assignee of the instant invention, and the disclosure therein is hereby incorporated by 
reference into the instant application. 

BACKGROUND OF THE INVENTION 

1 . Field of the Invention 

The present invention relates to a method for data communication of signals in 
units of a frame using an orthogonal frequency division modulation (OFDM) algorithm, 
and in particular, to methods of estimating carrier frequency offset error at the receiver. 

2. Description of the Prior Art 

OFDM system is a viable modulation scheme for data transmission time varying 
dynamic channels. However, it is known that performance of such system is highly 
susceptible to non-ideal synchronization parameters. Specifically, symbol timing and 
carrier frequency offset become an increasingly important issue in implementation of 
OFDM systems for practical applications. It is known that carrier frequency offset 
deteriorates performance of OFDM systems by introducing interference among the sub- 
channels. To overcome this imperfection, various compensation methods for estimation 
and correction of synchronization parameters have been proposed. In order to compare 
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the performance of these estimators, it is required to define a single number representing 
the goodness of the estimate. Assuming that all estimators are unbiased, i.e., expectation 
of the estimate is equal to the parameter, the variance of the estimator is used as a global 
measure for performance comparison of these estimators. 
5 Cramer-Rao lower bound (CRLB) is a fundamental lower bound on the variance 

of the estimators and the unbiased estimator whose variance equals CRLB is called 
efficient. When the evaluation of efficient estimator is not possible, it is desirable to 
obtain an estimator in which its performance becomes as close as possible to the CRLB 
fundamental bound. The estimator which is closest in performance to the CRLB 

10 estimator is known as a minimum variance unbiased (MVU) estimator. 

Categorically, the previously proposed methods for synchronization of OFDM 
systems can be classified into two main subclasses, namely minimum mean square error 
(MMSE) and maximum likelihood (ML) estimators. In MMSE approach, the estimator 
uses the information provided by the reference signal (pilot tones) in order to minimize a 

15 cost function associated with the synchronization parameters. A salient feature of this 
approach is that no probabilistic assumptions are made with regard to the data. Although 
MMSE estimators usually result in a tractable (globally stable) and easy to implement 
realization, no optimal criteria (probabilistic) is associated with these estimators. Also, 
since part of the transmitted information is allocated to the reference pilots, the 

20 bandwidth efficiency of these methods is lower in comparison to the non-pilot schemes. 

On the other hand, ML estimators provide the estimate of the unknown parameter 
subject to minimum probability of error criteria. Although not exactly efficient, ML 
estimators are asymptotically MVU, i.e., their variance attains that of MVU estimator as 
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the length of data record goes to infinity. However, due to the physical constraints, 
systems with infinitely long data records are not feasible for implementation purposes. 

P.H. Moose, in "A Technique for Orthogonal Frequency Division Multiplexing 
Frequency Offset Correction," in IEEE Trans. On Communications, Vol. 42, No. 10, pp. 
5 2908-2913, Oct. 1994, describes the use of a retransmission technique in order to reveal 
the frequency offset parameter in the likelihood function of the received signal. Due to 
the redundancy introduced by repeating the data block, the data rate efficiency is 
decreased by a factor of two. To avoid this imperfection, a ML estimator based on cyclic 
prefix (CP) is described by J. van de Beck, M. Sandel and P. O. Borjesson, in "ML 
10 Estimation of Timing and Frequency Offset in OFDM Systems," IEEE Trans. On Signal 
;fl Processing, Vol. 45, No. 3, pp. 1800-1805, July 1997. In this approach, the side 

■ 

information provided by the CP is used to obtain the likelihood function for joint 
estimation of symbol timing error and frequency offset in an OFDM system. 
]~ The likelihood function described in the Moose reference does not globally 

e E 15 characterize the observation vector over the entire range of the timing offset. 
H- Consequently, the ML estimator proposed based on this likelihood function would result 

O in considerable performance loss over a finite range of timing offset interval. 

Currently, there is increasing interest in multi-carrier modulation (MCM) for 
dividing a communication channel into several subchannels and transmitting many 
20 subcarriers through a single band using frequency division multiplexing (FDM) 

techniques. In the MCM method, however, because several subcarriers occupying a 
narrow frequency domain are transmitted at one time, a relatively longer symbol period 
results compared with a single carrier modulation method. The MCM method has, owing 
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to such characteristics, the advantages that equalization is easily performed and that it has 
immunity to impulse noise. OFDM is a type of the MCM designed to maximize the 
working frequency efficiency by securing orthogonality among the multiplexed 
subcarriers. OFDM is applied to mobile radio channels to attenuate multipath fading. 
5 In an OFDM transmitting/receiving system, modulation and demodulation of 

parallel data are carried out using the Fast Fourier Transform (FFT). It is required that the 
sampled data be sent in predetermined frames, having passed through a FFT routine, been 
time-division multiplexed, and transmitted, then restored at the receiving end. However, 
if the synchronization is in error in the course of restoring the frame, the signals 

10 demodulated after the FFT will be influenced by interchannel and intersymbol 
interference. Accordingly, the problem of synchronization in reforming the frame, 
especially any joint carrier frequency offset or symbol timing error, must be addressed as 
a matter of importance. 

Conventional synchronization methods as above-described encounter problems in 

15 that the process of synchronization is not only very complex, but the synchronization is 
not realized rapidly. 

SUMMARY OF THE INVENTION 
Motivated by the sub-optimum performance of estimators of the prior art, a 
likelihood function for joint estimation of carrier frequency offset and symbol timing 
20 error of OFDM systems is disclosed in the related, copending application; and a new 
optimum ML joint estimator is disclosed therein. In order to reduce the variance of that 
ML estimator, a new class of MVU estimators for frequency offset estimation of PFDM 
systems is disclosed herein. 
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There is disclosed to exist but one function of sufficient statistic which provides 
the MVU estimate of the frequency offset. The estimator provided by the instant 
invention is a closed form expression; providing an estimator which is a function of data 
statistic. Consequently, it does not suffer from converging to multiple local minima, a 
5 problem which arises in ML technique with nonconvex loglikelihood functions. 

The advantages of the instant MVU estimator over the class of previously 
proposed estimators are two; first, it is MVU, therefore its variance is minimum among 
the entire class of estimators which use the same probabilistic measure. Secondly, it 
provides a closed form expression for mapping the statistics into the estimation domain. 
10 The former property assures optimality of the estimator, while the later facilitate the 
closed loop analysis of the system. 

Accordingly, the present invention is directed at a synchronization method that 
substantially obviates on or more of the problems due to limitations and disadvantages of 
the prior art. To achieve these and other advantages, and in accordance with the purpose 
15 of the invention as embodied and broadly described, there is provided a method of 
estimating carrier frequency offset error in a received sample bit stream including an 
observation vector (OV), having an observed carrier frequency timing offset £, and a 
plurality of data-symbol frames, having a symbol timing offset error The method 
comprises the steps of generating a probability density function (PDF) based on the OV 
20 and generating from the PDF an estimate of carrier frequency offset error, £mvu I s> being 
a minimum variance unbiased (MVU) estimator. 

In another aspect of the invention, the OV comprises an L-bit cyclic extension 
portion and a first and a second N-bit synchronization frame, and wherein the PDF 
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comprises a first term, pi, based on the timing offset $ being within the span 1 to N and a 
second term, p2, based on said timing offset $ being within the span N+ 1 to N+L. In yet 
another aspect of the invention the received bit stream has uncorrelated independent 
identically distributed random signal and noise sequence variables with power of a] and 
&1 , respectively, wherein the OV is denoted x, and wherein the MVU estimator SmvuIs 
is the conditional expectation of a second moment estimator, said second moment 
estimator given by 
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In another aspect of the invention the MVU estimator smvu 1 9 is given by 
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where 3 is the imaginary operator and where 

S-N-\ IMA 



h=0 



1<3<N 
N+l<S<N+L 



The invention also provides a method of synchronizing a received sample bit 
stream, comprising the steps of transmitting at a transmitter the bit stream including an 
observation vector (OV), receiving and sampling at a receiver the bit stream, the sampled 
bit stream including, the OV with an observed carrier frequency offset e, and a plurality 
of data-symbol frames, having a symbol timing offset error generating a probability 
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density function (PDF) based on the OV; generating from the PDF an estimate of carrier 
frequency offset error, smvu 1 9, being a minimum variance unbiased (MVU) estimator; 
and synchronizing the received bit stream by the MVU estimate of carrier frequency 



Further aspects of the synchronization method of the invention include wherein 
the OV comprises an L-bit cyclic extension portion and a first and a second N-bit 
synchronization frame, and wherein the PDF comprises a first term, pi, based on the 
observed timing offset $ being within the span 1 to N and a second term, p2, based on the 
observed timing offset $ being within the span N+l to N+L. In yet another aspect of the 
inventuion, the synchronization method provides that the received bit stream has 
uncorrelated independent identically distributed random signal and noise sequence 

variables with power of a\ and a\ , respectively, wherein the OV is denoted x, and 
wherein the MVU estimator smvu I a is the conditional expectation of a second moment 
estimator, the second moment estimator given by 



Further aspect of the invention, the synchronization method wherein said MVU 
estimator smvu I s is given by 
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where 3 is the imaginary operator and where 
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\<3<N 
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In another aspect of the invention, A method of estimating carrier frequency 
offset error is provided for a received sample bit stream including an observation vector 
(OV) 5 having an observed carrier frequency timing offset e, and a plurality of data- 
symbol frames, the OV comprising an L-bit cyclic extension portion and a first and a 
second N-bit synchronization frame, comprising the steps of generating the expected 
value of the autocorrelation of the kth entry of the OV; and generating from the expected 
value an estimate of carrier frequency offset error, e m0 m> being a moment estimator. 

In yet another aspect of the invention, the moment carrier frequency offset error 
estimation method further includes wherein the OV is denoted x, and wherein the 



A method of synchronizing a received sample bit stream, is also provided by the 
invention comprising the steps of: transmitting at a transmitter said bit stream including 
an observation vector (OV); receiving and sampling at a receiver said bit stream, said 
sampled bit stream including, said OV with an observed carrier frequency offset e, and a 
plurality of data-symbol frames, having a symbol timing offset error $ ; 



moment estimator 




where the statistic T\(x) is defined as 
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generating the expected value of the autocorrelation of the Ath entry of the OV; 
generating from the expected value an estimate of carrier frequency offset error, e m om 5 
being a moment estimator; and synchronizing the received bit stream by the moment 
estimate of carrier frequency offset. 

Further aspect of the invention provides a synchronization method wherein the 
OV is denoted x, and wherein the moment estimator e mom is given by 



In yet another aspect of the invention, a method of estimating carrier frequency 
offset error in a received sample bit stream including an observation vector (OV), having 
an observed carrier frequency timing offset e, and a plurality of data-symbol frames, the 
OV comprising an L-bit cyclic extension portion and a first and a second N-bit 
synchronization frame, comprising the steps of: selecting one estimator method from a 
group of three: a maximum likelihood (ML) estimator, a minimum variance unbiased 
(MVU) estimator, a moment estimator, based on said OV; and generating from the 
selected estimator an estimate of carrier frequency offset error. 

Additional features and advantages of the invention will be set forth in the 
description which follows, and in part will be apparent from the description, or may be 
learned by practice of the invention. The objectives and other advantages of the invention 
will be realized and attained by the methods particularly pointed out in the written 
description and claims hereof, as well as the appended drawings. 




where the statistic Tj(x) is defined as 




BRIEF DESCRIPTION OF THE DRAWINGS 
FIG. 1 is a signal flow block diagram depicting the frequency recovery loop for 
estimation of frequency offset. 

FIG. 2 is a block diagram depiction of the ML estimator disclosed in the cross- 
referenced, copending application and the MVU and moment estimators of the instant 
invention. 

FIG. 3 is a graph comparing the computer simulations of the prior art ML 
frequency offset estimation and the ML frequency offset estimation disclosed in the 
related, copending application. 

FIG. 4 is a graph comparing the computer simulations of the MVU frequency 
offset estimator of the instant invention with that of the CRLB estimator. 

FIG. 5 graphically depicts the closed loop performance of the MVU estimator of 
the instant invention by comparing the frequency offset estimate of the closed loop MVU 
estimator against an analytical derivation thereof. 

DETAILED DESCRIPTION OF THE INVENTION 

The cross-referenced co-pending application entitled "Globally Optimum 

Maximum Likelihood Estimation of Joint Carrier Frequency Offset and Symbol Timing 

Error," U.S. Ser. No. 09/496,890, filed on March 2, 2000, incorporated herein by 

reference, discloses a probability density function (PDF),p, which globally characterizes 

an observation vector x according to the equation 

p(x y e 9 9) = p, (x 9 e 9 9)(U[9 - 1]- u[S - * - 1]) (U) 
+ p 2 (x i e , 3)(U[& - N - 1]- U[&- N ~L + 1]) 
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where e and S are the frequency offset and symbol timing error introduced by the 
synchronization mismatch in the carrier frequency and symbol timing, respectively. And 
pi and P2 are probability density functions derived and described in the related copending 
application. 

5 Further derived and described therein are the statistics Ti and T2 



and 



4V-1 N+1A 



(15) 



k=0 k£ 



T 2 (x,0) = < 



£ .|jc[jfc]| 2 + |x[jt + tf]| 2 1<9<N 

i-AT-1 , AT+£-l 

£|jc[*]| +|jc[* + AT]| + £ +|jc[fe + iV]| N + \<3<N + L 



(16) 



10 N and Z being the number of samples per frame and the number of cyclically-extended 
samples appended at the beginning of each frame, respectively, all as further described in 
the related, copending application. 

The MVU estimator of the instant invention is derived by resorting to the theory 
of sufficient statistics. The first step in deriving the MVU estimator is to obtain the 

15 sufficient statistic for the PDF given in (14). The sufficient statistic is known to be a 
function of the observation vector, namely T(x), such that the conditional PDF of the 
observation vector given T(x) does not depend on the unknown estimation parameters 
[€,&]. Evaluating the sufficient statistic is a formidable task for the broad class of PDFs, 
however the Neyman-Fisher Factorization theorem can be used for identifying the 

20 potential sufficient statistics. According to this theorem, if the PDF can be factored in the 
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form g(T(x\€,9,)h(x) where g is a function depending on x only through T(\) and h(x) 
is a function depending only on x, then T(x) is a sufficient statistic for estimation of the 
parameters s and 9 . By reformulating the PDF given in 14 to 

24e^rJx^)|a|]Ha| 2 r 2 (x jl 9) 



2^1-|a| 2 ^crr +cr 



2 2 

5 « 



p(x,e,9)=e \ A s n ) /^( x ) 
there is a direct dependency between the parameter 9 and the statistics Ti( x , 9 ) and 
r 2 (x ? 9) . Based on this observation, the Neyman-Fisher theorem fails to provide a 
sufficient statistic for estimation of 9 , we can factor the PDF into 

dy^^M) |a| 



p(x 5 ^| i9 )=e ^ a ^ hl {x) (17) 

Clearly then, ^(x,^) forms a sufficient statistic for estimation of the parameter s . 

Next, application of the Rao-Blackwell-Lehmann-Scheffe (RBLS) theorem is 
used to find the MVU estimator. According to this theorem, if e is an unbiased 

estimator of s and r(x) is a sufficient complete statistic for € then s = E(e | r(x) ) is a 

valid, unbiased, MVU estimator of s . 

In applying the above theorem, we need to obtain an unbiased estimator of 
s, termed s , and determine the conditional expectation of this estimator given the 
statistic T x (x, 9-). An appropriate candidate for the unbiased estimator of £ can be 

obtained from the statistical moments of the random vector x. The autocorrelation 
function of the observation vector x, E, yields that the second moment of the random 

variable x[k] with k € Q, satisfies the following identity 
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E^c[k\c*[k + N\}=(7h~ i27tS (18) 
Having this observation, one uses the second moment estimator as an unbiased estimator 
for s as given by 



s = —\n\ 
in 



1 « 9 + Z " 1 r i 

1 £ Jdjfc]x*[£ + ;V] 



W k = 9 



(19) 



where 3 is the imaginary operator. In deriving the above estimator, *[k + JV]) 

was replaced by its natural estimator — V 9+1 ' x[k]x * [k + n] It is straightforward to verify 

J. *k--3 



that this estimator is unbiased as it satisfies the condition 

•9+1-1 



E{e} = — InE 
2n 
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(20) 



Next, one obtains the conditional expectation of € given the sufficient statistic as 
follows 



€ MVUl^ 
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(21) 



It is important to emphasize that since the underlying PDF given in (14) belongs 
to the exponential family of PDFs, then the sufficient statistics ?i(x,i9) forms a complete 
statistic for estimation of the parameter s . Therefore, the mapping function obtained 
from applying RBLS theorem, namely ln7(x, <9), is but one function of the statistic 
T x (x, &) and no other estimator with the same statistic can result in a lower variance with 
respect to MVU estimator. 
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A. Cramer-Rao Lower Bound 

It is known that under broad conditions, the variance of any unbiased estimator of 
a nonrandom parameter s satisfies the CRLB as 



var 



(22) 



where l{s) is the Fisher Information given by 

d 2 \npr(*,e\ 3 ) 



ds 2 



(23) 



Substituting (17) into (23), after some algebraic manipulations, the CRLB of the MVU 
estimator becomes 



1 + 



1 



SNR 



-1 



{ln) 2 L 



(24) 



10 where SNR ^- is the signal to noise ratio at the receiving end. 

B. Closed Loop Performance 

FIG. 1 is a signal flow block diagram depicting the frequency recovery loop for 
estimation of frequency offset according to the instant invention. A closed loop system is 
obtained by feeding back the information obtained from the estimator into a sampler 
15 block (bootstrap). A sampler 10 updates its frequency at the beginning of each 

observation vector (every (N+L) samples). The sampler 10 receives the latest signal and 
produces the sampled symbol stream x[k] therefrom which is conducted to a block 12 and 
a block 14 performing the autocorrelation operation on the observation vector as required 
by equation (19). The result of the autocorrelation operation is conducted to a moving 
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average (MA) filter block 16 and therefrom to a block 18 performing the natural 
logarithm function and thence to a block 20 performing the imaginary operation function, 
as required by equation (19). To match the various sampling frequencies used in the 
system, a down sampler (decimator) block 22 is used prior to return to the sampler 10 to 
produce the minimum variance unbiased estimator, s [m]. Finally, a gain block (G) 24 is 
used to control the closed loop characteristic of the system (stability, settling time, noise 
sensitivity). 

According to Fig. 1, the frequency offset for the mth observation vector can be 
expressed as 

where Az[m - \]={£[m - \]-e)[m - 2])G and x m + z)+ 5 + /} The term inside the sum is 

a stochastic quantity and does not have a closed form expression. However, for 
reasonably high signal to noise ratio it can be well approximated by its expected value 
15 H*^,.,) [ / l x * (m _2) M] = al s ) • Th ere f" ore > the expression inside the sum can be written as 

L - X . 2Mm~m , -J )r J - f 2*fr»-lfc-l) SU1 — 



/=0 i ^ J n sin — 1 J 



IN ) 

_ 2nLe\m-\\L-\) L 

As[m-l]-»0 e J 2N (26) 
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Substituting (26) into (25), after some algebraic manipulations, the frequency offset of 
m'th observation vector becomes 

£H=-(5[ w -l]-4n-2])^(5+^) (27) 

The above equation represents a second order finite difference system in which its 
dynamic can be obtained from solving the following equation 
s[m] + ps[m - 1] - /3s[m - 1] = 0 

where p= ^3 + ~y~ j • Clearly, the solution to the above finite difference equation is the 

form of 

e[m) = c x {z x )>» + c 2 {z 2 f (28) 



-p±4p 2+ ap 

where Zj 2 = 2 are two dynamical modes of the system. The smaller root 

(negative) results in a high frequency oscillation in the frequency offset estimate. 
However, as is shown in a computer simulation, infra , this term is filtered out by the 
moving average filter. To assure stability, the gain block should be set such that both 
poles lie inside the unit circle. 



■p m + V/? 2 max+4/?max|< 2 (29) 



where 



A G 
/?„,,= max — 
/>max 9 N 



G_ 
N 



N + L + — I (30) 



2 ; 

It is known that when the timing offset parameter is not known to the receiver or 
if the noise PDF differs from Gaussian distribution, finding the optimum estimator (ML, 
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MVU, CRLB) may not be an easy task. However there exists a moment estimator which 
provides a consistent estimate for estimation of frequency offset regardless of noise 
distribution and timing offset values. Although there is no optimum criterion associated 
with the moment estimator due to its simple structure, it is frequently used as an initial 
estimate for other estimators such as ML estimator. Consider a sequence of first N+L 
samples of vector x. Using the autocorrelation of klh entry of this vector satisfies the 
following identity Moment Estimator 



rxxW = 



2 - jlns keQ 
s 

o ken 



(31) 



Using Base Rule, the expected value of the above function (with respect to parameter k) 
can be expressed as 

(32) 



E k 
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xx 



[N] 



= cr} exp ">2» pr (ke£l)+0pr (k Z Cl )= -r^rr-a } exp -J*** 

l* "\~ yv 

Substituting the Mh autocorrelation lag with its natural estimator, the moment estimator 
for frequency offset under uncertain timing offset can be found as 

imom=^T 3 (x)} (33) 

where the statistic T^(x) is defined as 

JV+Z+1 (34) 



3 Lo? t=0 J 



k=Q 

Statistical assessment of moment estimator is a formidable task over entire range of SNR. 
However, for relatively high SNR, the random observation vector is heavily concentrated 
about its mean. Using the statistical linearization, and a first-order Taylor expression of 
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the estimator about its mean to obtain the variance of estimate. In doing so, one 
substitutes for the random variable x in (34) and obtains 

Inks 



Ti([x])=/(M)= f{^5 2 "l< N 



k=0 



2n{k+N]e r ,1 
s*[k-3+N]exp N L J 



}(35) 



where the signal (s) and noise (w) vector are defined as 

s A [s[6^ --s[n + L-i\] 
wAMO]]...M{iV + Z-l]] 



(36,37) 



By virtue of the above equation, the expected value of observation vector for a 
fixed realization of signal vector s would be E\s mom ]= f(s) . Then performing a first 
order Taylor expansion of f([s.w]) about the point £>v[ x ] yields 

smom =/([^0])+VM/([5,0])*>v=/(5,0]) + Z ' f l" 1 ^%^{w = 0)^ (38) 



Taking the derivative of (35) with respect to wfnj and setting w=0, results in 



dh 



dM[n] 



1 N+L-l 



w=0 = 



( S s[i-3]s*[i-3 + N]e j2n€ ) _1 s*[n + N-d]e N 



2tt /=0 



The second term in (38) represents the contribution of noise in the estimate. 

Knowing that noise samples are iid with power of a\ , the variance of estimate can be 
obtained from 



\3r(s ) = 



2 2 



n=0 



< N+L-\ 



(39) 
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For sufficiently large block lengths (AO, the above term can be well approximated as 



The resemblance between the estimators of equations (21) and (33) reveals a 
unified structure which characterizes the ML, MVU and moment estimators which can be 
classified into a single unified structure. This provides a unique framework for analysis 
of the proposed estimators. Moreover, it allows an investigation of the effect of symbol 
timing error in the estimation of carrier offset for each individual estimator. Comparing 
the MVU estimator given in (21) to the moment estimator in (33) reveals some 
similarities in the structure of the estimators. Clearly, both moment and MVU estimators 
use the same mapping function, namely the log function, to project the data statistics into 
the estimation domain. The only difference is in the form of statistics used for each 
scheme. 

FIG. 2 is a block diagram depiction of the ML estimator disclosed in the cross- 
referenced, copending application and the MVU and moment estimators of the instant 
invention. An autocorrelation operation on the observation vector x[k] is performed by 
blocks 30 and 34; as shown in FIG. 2, both estimators of the instant invention obtain the 
statistic by correlating the samples with the Mh delayed samples. This operation is 
performed by using a moving average (MA) filter (shown as a block 36) in the structure 
of estimators. However, the MVU and moment estimators use different upper and lower 
bound for the MA filter. In the moment estimator, the averaging is performed over the 
first N+L samples of the observation vector. This would remove the requirement of 
knowing the exact timing offset parameter in estimation of carrier frequency offset. 




(N + L)cr}al_ (N + L) 



(40) 




19 





LTI014 

However, the estimate obtained from using this estimator results in a less accurate 
estimate (more variance) in comparison to MVU estimate. On the other hand, the MVU 
estimator requires the knowledge of symbol timing in the estimation of carrier frequency 
offset, as shown in FIG. 2 by a block 38 conducting S to the MA block 36. Finally, as in 
5 FIG. 1 , a block 40 and a block 42 perform the natural logarithm and imaginary, 

respectively, operations on the symbol stream before generating the estimate of carrier 
frequency offset 8. 

Although the resemblance between MVU and ML estimators may not be as 
evident as that of MVU and moment estimator, it can be shown that ML estimator can 
10 also be classified into the same family. Knowing the fact that sjlogjT^ ^j}= ZTfaS) the 

ML estimator can be expressed as 



*u, ^4 7 s* + [*+4 (4i) 

ML In [k=a J 



where the parameters a and /?are functions of 3 ml and can be obtained from the 
expression for Ti(x,d), supra . Thus, the ML estimator falls into the same family of 

15 estimators. FIG. 2 also displays in tabular form these estimators, their respective 
summation upper and lower bounds, and their respective averaging intervals. It is 
noticed that the ML estimator provides the upper and lower bound of the moving average 
filter (36) by extracting the timing parameter from the likelihood function. Although the 
ML estimator has the advantage of exploiting the entire bandwidth by removing the 

20 requirement for having pilot tones, the symbol timing estimate obtained from the ML 

estimator has a larger confidence interval. This may result in a considerable performance 
degradation in comparison to the pilot-based schemes. 
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Computer simulation is used to assess the performance of the three estimators for 
synchronization of an OFDM system. The variance of estimator is used as a performance 
measure through the study. The simulation parameters used are typical of the OFDM and 
digital audio broadcast (DAB) environments. More specifically, the chosen FFT size (AO 
5 for OFDM is 64. Unless otherwise specified, the length of cyclic prefix (L), signal to 
noise ratio, and frequency offset are set to 8, 20 dB and .01, respectively. Monte Carlo 
simulation is used to evaluate the performance of the three estimators. 

FIG. 3 graphically depicts a comparison between the performance of the ML 
frequency offset estimator disclosed in the related, copending application with the prior 
10 art ML frequency offset estimator given in the van de Beek, et al, reference over the 
1% range of timing offset parameter (& e [l, N + L] ) 

FIG. 4 graphically depicts the performance of the MVU frequency offset 
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\n estimator of the instant invention under complete knowledge of timing offset error. A 
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; ^ careful examination of the variances reveals that the gap between MVU estimator and 

!*! 15 CRLB tends to zero as SNR increases. Also as illustrated in FIG. 4, the departure from 
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;Z CRLB happens rapidly as SNR goes below a threshold. The threshold also depends on 
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w the length of CP and is moved toward lower SNRs as L increases. This can be justified in 

terms of having more observation samples in estimating the unknown parameter. The 
choice of cyclic prefix length L represents a tradeoff between data rate reduction and 
20 performance (lower variance). Increasing L brings the performance of MVU estimator 
closer to the CRLB, nevertheless, it could result in a considerable data rate reduction due 
to the redundancy introduced by CP. 
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